%calculates the local time from fday and longitude
%also calculates the LT at India from longitude
%latest date 31 OCT 2005 (B'day !)

function[lt1,lt2] = cham_lt(fday,lon),

if lon < 0,
    lon = lon + 360;
end;


if fday < 0,
    lt1 = (fday-floor(fday) + 1) * 24 ;
else,
    lt1 = (fday-floor(fday))* 24;
end;
    
    lt1 = lt1 + lon/15;
    if lt1 > 24,
        lt1 = lt1-24;
    end;
lt2 = (fday-floor(fday))*24 + 5.5; % 78 is the longitude of ETT

    if lt2 > 24,
        lt2 = lt2-24;
    end;
